Study of the pvr Puzzle in Charmonium Decays 



MO Xiao-HuB YUAN Chang-ZhengQ and WANG Pin^ 

Institute of High Energy Physics, CAS, Beijing 100039, China 
(Dated: February 2, 2008) 

The theoretical explanations about the "pvr puzzle" in charmonium decays are reviewed exten- 
sively, and the comparison of theoretical predications with experimental data is made whenever 
possible. Three methods to estimate the ratio of the branching fractions of J /tj) and t/i' decays are 
also discussed. It is pointed out that in order to understand the pn puzzle, and the dynamics of 
charmonium decays, systematic studies should be made in theory, phenomenology and experiment 
aspects. 

PACS numbers: 13.25.Gv, 12.38. Qk, 14.40.Gx 



I. INTRODUCTION 

Crisply defined experimental puzzles in high-energy 
physics always have the prospect of leading to new discov- 
eries. One prominent example is the 9-t puzzle of 1956 
which led to the parity revolution. Therefore puzzles in 
physics often draw considerable attention to theorists. 
The ratios of hadronic decays of the ?/'(3686) (shortened 
as ij}') to the same decays of the J /ip is such a puzzle 
which has been studied substantially since 1983. 

Since the OZI suppressed decays of J/-0 and V' to 
hadrons are via three gluons or a photon, in either case, 
the perturbative QCD (pQCD) provides a relation [l| 

Qh = #^ = « 12.7% . (1) 

This relation is referred to as the "12% rule" which is 
expected to be held to a reasonable good degree for both 
inclusive and exclusive decays. The so-called "/ott puzzle" 
is that the prediction by Eq. ^ is severely violated in 
the piT and several other decay channels. The first evi- 
dence for this effect was found by Mark-II Collaboration 
in 1983 0. From then on many theoretical explanations 
have been put forth to decipher this puzzle. 

With the recent experiment results from BESII 
and CLEOc about J/^ and -0' two-body decays, 
such as vector-pseudoscalar {VP), vector-tensor {VT), 
pseudoscalar-pseudoscalar (PP), and baryon-antibaryon 
{BB) modes, and about multi-body decays at the JM, 
the ip' or even at the i/'(3770) (shortened as ip") a 
variety of solutions proposed for the puzzle can be tested 
at the level of higher accuracy. In this treatise, we survey 
the theoretical works on the pir puzzle and compare them 
with the available experimental data. From the theoret- 
ical point of view, since the Q-value is smaller than 12% 
for pTT, it may be caused either by enhanced or suppressed 
J/ip decay rate. Another possibility is by both. So we 
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classify the relevant theoretical speculations into three 
categories: 

1. J/i/j-enhancement hypothesis, which attributes the 
small Q-yalue to the enhanced branching fraction 
of J/tjj decays. 

2. ?/;'-suppress hypothesis, which attributes the small 
Q-value to the suppressed branching fraction of ijj' 
decays. 

3. Other hypotheses, which are not included in the 
above two categories. 

In the following content, first reviewed are the theoreti- 
cal works on pn puzzle, and the predictions from them are 
compared with the newly available experimental results; 
then expounded are three methods to estimate the ratio 
between -0' and J/t/j decays into the same final states; 
after that some comments are made on the implications 
from the review on pir puzzle; at last there is a short 
summary. 



II. REVIEW OF THEORETICAL WORKS ON pn 
PUZZLE 

A. J/?/)-enhancement Theory 

In the earlier days of the pn puzzle, it was noticed 
that the decay of 1 charmonium into pir final state vio- 
lates the Hadronic Helicity Conservation (HHC) theorem 
(see below for expound) |24i], and so such decay should 
be suppressed. Therefore people think there must be 
some mechanism which leads to the great enhancement 
for J/ip — > pn decays. The two schemes presented in this 
section were proposed following this line of reasoning. 

1. J /tj)- glueball Admixture Scheme 

The idea of J /'ip decays via a glueball was proposed 
by Freund and Nambu [2^ (FN hereafter) soon after 
the discovery of J particle to explain the breaking 
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of Okubo-Zweig-Iizuka (OZI) rule [2y|. In such mech- 
anism, the breaking results from the mixing of the w, 
0, and J /i}} mesons with an SU (4)-singlet vector meson 
O. They found that such an O meson should lie in the 
1.4-1.8 GeV/c^ mass range with the width greater than 
40 MeV/c^, and it should decay copiously into p7r, K*K 
while exhibiting severe suppression of decays into KK, 
e^e~ and modes. These authors presented several 

quantitative predications for experimental search. Two 
of them are 



Ri = 



= 0.0115-0.087 



i?2 = 



JH^^ < 8 X 10-5 . 



With the current available data and using the three pi- 
ous final state as a substitute for p7r in both cf) [23| and 
J/V' ll, Hi decays, we obtain the first ratio i?i ~ 0.003, 
which is almost one order of magnitude smaller than 
the predication. For the second ratio, by virtue of 
PDG [271 value for K~^K~ and new experiment result for 
K^Kl [6], it is estimated that B(J/ij KK) - 10"^^, 
together with the results for p-K we have i?2 ~ 10^^ 

which is much larger than the predication. 

The first attempt to explain the pir puzzle in terms of a 
glueball near J/ip was proposed by Hou and Soni [28i] (HS 
hereafter). They attributed the enhancement oi J/ip 
K* K and J /ijj pvr decay modes to the mixing of the ij: 
with a J^'-' — vector gluonium, also designated by 
O. The differences between FN's and HS's pictures lie in 
the following aspects: 

• Based on the potential model applied to glueball, 
the mass of a low-lying three-glue state is estimated 
to be around 2.4 GeV/c^ [1^, rather than 1.4 to 1.8 
GeV/c^ in Ref. [H. 

• The mixing of O with ip' is taken into account, 
which has been ignored in previous work. 

• Since the gauge coupling constant in QCD is mo- 
mentum dependent, the mixing parameter is taken 
to be a function of the invariant mass , which de- 
creases rather sharply with the increase in . Such 
propagator effect gives rise to a more suppression 
on the decay rates of ip' relative to J jij] for decays 
to p-K and K*K channels. 

By virtue of their assumption, HS suggested a search 
for the vector gluonium state in certain hadronic decays 
of the -0', such as V'' — > tttt -I- X, -t- X, ?/ -f X , where X 
decays into VP final states [28] . 

Based on HS's idea, Brodsky, Lepage, and Tuan [s^l 
(BLT) refined the glueball hypothesis for the pi: puzzle. 
They assumed the general validity of the pQCD theo- 
rem that the total hadron helicity is conserved in high- 
momentum-transfer exclusive processes, in which case 



the decays to pir and K*K are forbidden for both the J/V' 
and V''- This pQCD theorem is often referred to as the 
rule of Hadronic Helicity Conservation (HHC) [13], which 
is based on the assumption of the short-range "pointlike" 
interactions among the constituent quarks throughout. 
For instance, J/'ii){cc) — > 3g has a short range ~ 1/toc 
associated with the short time scale of interaction. Nev- 
ertheless, if subsequently the three gluons were to res- 
onate forming an intermediate gluonium state O which 
has a large transverse size covering an extended time pe- 
riod, then HHC would become invalid. In essence the HS 
model takes over in this latter stage. 

Final states h which proceed only through the inter- 
mediate gluonium state satisfy the ratio 



Qh = 



B{'4'' ^ e+e-) (Mj/^ - Mp)^ + ry4 
B{J/^P ^ e+e-) (M^, - MqY + r2,/4 



(2) 



The Qh is small if the mass of O is close to the mass 
of J/V'- The experimental limits at that time [1, HO, HH 
imply that the O mass is within 80 MeV/c^ of the mass of 
J /if: and its total width is less than 160 MeV/c^. Brodsky 
et al. recommended a direct way to search for O, that 
is to scan the e+e" VP cross section across the J/ip 
resonance. 

Another related work by Chan and Hou [s^l studied 
the mixing angle 0oii> and the mixing amplitude /o^ of 
the J/tp and vector glueball O based on the framework 
of potential models of heavy quarks and constituent glu- 
ons. They obtain |tan6'c)^| = 0.015 and /o^a(™c)i/>) ~ 
0.008 GeV^ 

On the experimental part, BES has searched for this 
hypothetical particle in a pi: scan across the J/4' mass 
region in e+e^ annihilations as well as in the decays of 
-0' —^ ttttO, O — > pTT, and found no evidence for its ex- 
istence [1, 111] . The data constrains the mass and width 
of the O to the range | Mo - Mj/^ \< 80 MeV/c^ and 
4 < To < 50 MeV/c^ Although the absence of dis- 
tortion in BES energy scan of J/0 pir does not rule 
out Mo ~ Mj/^, it puts a lower bound to To- However, 
as indicated in Ref. [3§], the experimentally constrained 
mass is several hundred MeV/c^ lower than the mass of 
the lightest vector glueball calculated in lattice simula- 
tions of QCD without dynamical quarks [STj . 

Recently, the experimental data from BES and CLEOc 
turned out to be unfavorable to this glueball hypothesis. 
Among them is the observed large branching fractions 
of the isospin- violating VP mode -0' cun^ [l3i EBl ■ 
This contradicts with the assertion that the pattern of 
suppression depends on the spin-parity of the final state 
mesons. In addition, according to BLT's analysis, one 
obtains the relation l38l 



B{J/^p ^ Lun°) 
B{J/iIj p"7rO) 



< 0.0037 



which is much smaller than the PDG06 value 0.08 [27| . 
Another experimental result which is unfavorable to such 
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hypothesis is the suppression of tp' decays into vector- 
tensor (VT) final states Since hadronic VT de- 
cays, unlike the VP decays, conserve HHC, some other 
mechanisms must be responsible for this suppression in 
the model. Furthermore, it has been argued that the O 
may also explain the decay of J/ip into 0/o (named pre- 
viously S*) but not to /9ao(980) (named previously 6), 
since the O mixes with the (j) and enhances a mode that 
would be otherwise suppressed [s^l- However, the decay 
tp' — > (j)fo [13] is not suppressed as found by experiments, 
which implies the absence of anomalous enhancement 
in J ftp — > (pfo, thus contradicts with this explanation. 
Anselmino et al. extended the idea of J/ip-0 mixing to 
the case of rjc VV and pp [3§| . They suggested that 
the enhancement of these decays can be attributed to the 
presence of a trigluonium pseudoscalar state with a mass 
close to the rjc mass. So far there is no experimental 
evidence for such a state. 

In fact, as pointed out in Ref. [13], this glueball ex- 
planation has some unanswered questions: (i) Why only 
the pir and K*K channels are affected but not O — *■ Stt 
etc. viz., one must assume in an ad hoc way that the O 
couples predominantly to pir and K*K; and (ii) if such 
a narrow heavy 1 gluonium state exists, why have the 
narrow 2++ states not been found, which should be 
lighter and easier to search? 

2. Intrinsic-charm-component Scheme 

Brodsky and Karliner (BK) put forth an explana- 
tion for the puzzle based on the existence of the in- 
trinsic charm \qqcc) Fock components of the light vector 
mesons [4l| . They noticed the fact that quantum fluctua- 
tions in a QCD bound state wave function will inevitably 
produce Fock states containing heavy quark pairs. The 
intrinsic heavy quark pairs are multiconnected to the va- 
lence quarks of the light hadrons, and the wave functions 
describing these configurations will have maximal ampli- 
tude at minimal off-shellness and minimal invariant mass. 
In the case of the p meson, with the consideration of the 
light-cone Fock representation: 

p+ = ip'^-\ud) + ip''- _\udcc) -I- • ■ • . 

ud^ ' udcc^ ' 

Here we expect the wave function of the cc quarks to be 
in an S'-wave configuration with no nodes in its radial de- 
pendence, in order to minimize the kinetic energy of the 
charm quarks and thus also minimize the total invariant 
mass. 

The presence of the \udcc) Fock state in the p allows 
the J /%p pi: decay to proceed through rearrangement 
of the incoming and outgoing quark lines; in fact, the 
\udcc) Fock state wave function has a good overlap with 
the radial and spin |cc) and \ud) wave functions of the 
J ftp and pion. On the other hand, the overlap with the 
Tp' will be suppressed, since the radial wave function of 
the n = 2 quarkonium state is orthogonal to the nodeless 
cc in the \udcc) state of the pir. Similarly, the \uscc) Fock 



component of the K* favors the J/tpK configuration, al- 
lowing the J/ip K*K decay to also proceed by quark 
line rearrangement, rather than cc annihilation. 

These authors also suggested comparing branching 
fractions for the r]c and 77^ decays as clues to the im- 
portance of rjc intrinsic charm excitations in the wave- 
functions of light hadrons. 

B. i/)'-suppress Theory 

The hypothesis of the existence of a glueball to explain 
the pTT has been questioned soon after it was proposed. 
In addition, it is also pointed out [i^l that the helicity 
suppression is not a strong constraint in the charmonium 
energy scale. Under such case, one comes naturally to 
the idea that it is not J /ip — > p-n which is enhanced, but 
rather ip' — > pi: which is suppressed. Seven explanations 
or models collected in this section are put forth along this 
line. 



1. Sequential- fragmentation Model 

Karl and Roberts have suggested explain the pn puz- 
zle based on the mechanism of sequential quark pair cre- 
ation [4^. The idea is that the quark-antiquark pairs 
are produced sequentially, as a result the amplitude to 
produce two mesons in their ground state is an oscilla- 
tory function of the total energy of the system. They ar- 
gue that the oscillatory fragmentation probability could 
have a minimum near the mass of ip' , which provides an 
explanation for the suppressed ip' decay. Even though 
their evaluations could generally accommodate the data 
for decays of J /ip and ip' to pir and K*K, it runs into 
difficulties when it is extrapolated to T decays. Accord- 
ing to their calculation, the oscillations of probability 
amplitude are damped out in the region of the T res- 
onances, so the pn channel is present in the decay of all 
T, T', T", • • • resonances with a common rate. This 
leads to a prediction r(T pn) = 0.05 keV, or equiv- 
alently B{T pir) = 9.4 x 10^^, which is above the 
current upper hmit B{T ^ pir) < 2 x 10"" [13]. More- 
over, their calculation seems hard to explain the large 
branching fraction for (p decays to pn 27] due to the fact 
that their fragmentation probability tends to zero as the 
mass of the piv decayin g sy stem approaches IGeV/c^. 

In a further analysis [4J|, Karl and Tuan pointed out 
that if a suppression is observed in three-meson channels 
the explanation based on sequential pair creation would 
be undermined. Recently such a suppressed channel, viz. 
(PKK, is found by CLEOc [2l|. 

2. Exponential-form-factor Model 

Guided by suppressed ratios of ip' to J /ip decays to 
two-body hadronic modes, Chaichian and Tornqvist sug- 
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gested [iol] that the hadronic form factors fall exponen- 
tially as described by the overlap of wave functions within 
a nonrelativistic quark model. This behavior explains 
the drastically suppressed two-body decay rates of the 
ip' compared with those of the J/tp. Recently, the re- 
port on observation of a number of VP channels in tp' 
decays 0, M, El 

such as LOT]' , (prj' , pr( has proved that 
the predicted decay fractions are overestimated. More- 
over, the branching fraction for u^tt" [23|, is well below 
the prediction by this model which is 1.04 x 10^^. 

Another problem of the model is that it does not single 
out just the VP channel, the other channels, for example 
VT channel, are also estimated to have small branching 
fractions which are not compatible with the BES mea- 
sured results ligl. 



3. Generalized Hindered Ml Transition Model 

A so-called generalized hindered Ml transition model 
is proposed by Pinsky as a solution for the puzzle [45| . 
It is argued that because J /ip — > 777 is an allowed Ml 
transition while ip' —> jrj' is hindered (in the nonrela- 
tivistic limit), using the vector-dominance model to relate 
Ip' — > 777' to ip' ipv' one could find the coupling G^i^r^^ 
is much smaller than G^-^rj^ , and then by analogy, the 
coupling Guj'p-K would be much smaller than G^jpTi- Here 
GujpTT can be extracted from data by virtue of the anal- 
ysis using the vector-dominance model and a standard 
parameterization of OZI process (46j . Then assuming 
■0' —>■ /CTT to proceed via ip'-uj' mixing, while J/ip — * p7r 
via ip-LO mixing, one would find that ■0' — > pir is much 
more severely suppressed than J/tp — > pn. The similar 
estimation could be preformed for K*K and other VP 
final state, and one can expect a suppressed Q: 



B{iP' VP) 
B{tP 



= 1.47 



TtoM) (G 



TV'VP 



Fv 



0.06% 



VP) ■ Ttotii'') V Gvvp 

(3) 

where Fyi I Fy = 0.3, G p-n I G p^^ — 0.066 according 
to Ref. [43 • This Q is much smaller than the present 
experimental results [1, EB] . 

Moreover, in this model, the coupling G^'^f^ for w' 
ujf2 should not be suppressed because by analogy the 
couphng G'^'^xc2 is not small due to the fact that the El 
transition ip' 7Xc2 is not hindered Therefore via 
tp'-uj' mixing the ip' uj' —f Lof2 decay is expected to be 
not suppressed, which contradicts the BES result [l9|. 



4- Higher- Fock- state Scheme 

Chen and Braaten (CB) proposed an explanation [H 
for the jOTT puzzle, arguing that the decay J /ip — > pi: is 
dominated by a Fock state in which the cc pair is in a 
color-octet ^5*1 state which decays via cc — > qq, while the 
suppression of this decay mode for the ip' is attributed to 
a dynamical effect due to the small energy gap between 



the mass of the ip' and the DD threshold. Using the 
BES data on the branching fractions into pn and K*K 
as input, they predicted the branching fractions for many 
other VP decay modes of the ip' , as listed in Table |T1 from 
which we see most measured values falling in the scopes 
of predictions, but we also note for cjtt mode, the devi- 
ation from the prediction is obvious. Here it should be 
noticed that the values adduced in Table U are calculated 
on the strength of the measured branching fractions from 
earlier experiment, the new measurements on the branch- 
ing fractions for pn and K*°Tf + cc. from BES [ll|. El 
and CLEOc [ll may have impact on the corresponding 
evaluations. 



TABLE I: Predictions and measurements for Qvp in unit of 
1% for all VP final states. The value for p-K and K*^K -\- cc. 
from Ref. [isl ] were used as input. The theoretical parameter 
X — 0.64 is due to results of Ref . 14911 and the experimental 
results come from Ref. [lol . [Til . |5C|. 



VP 


X — 


: 0.64 


Exp. 


pvr 


- 


0.25 


0.13 ±0.03 


K*'^JiO + c.c 


1.2 


-3.0 


3.2 ±0.08 


K*+K- +C.C 


0- 


0.36 


"^iQ+o ^^ 

^•0^-0.36 


UJTj 


0- 


- 1.6 


< 2.0 


UT)' 


12 


- 55 


J^^-13 


<Pv 
H 


0.4 


-3.0 


5.1 ± 1.9 


0.5 


- 2.2 


9.4 ±4.8 


PV 


14 


- 22 


g 2 + 3.6 


pri' 


12 


- 20 


17'.8+lft 




11 


- 17 


44+1.9 



Besides the predictions in Table HI CB's proposal also 
has implications for the angular distributions for two- 
body decay modes. In general, the angular distribution 
has the form 1 + acos^^, with — 1 < a < ±1. CB's 
conclusion implies that the parameter a for any two-body 
decay of the ip' should be less than or equal to that for 
the corresponding J/ip decays. But this needs further 
supports from analyses based on large data sample in 
the future. 



5. Survival- chamonia- amplitude Explanation 

A model put forward by Gerard and Weyers enter- 
tains the assumption that the three-gluon annihilation 
amplitude and the QED amplitude add incoherently in 
all channels in J/ip decays into light hadrons, while in the 
case of Ip' decays the dominant QCD annihilation ampli- 
tude is not into three gluons, but into a specific config- 
uration of five gluons [5l| . More precisely, they suggest 
that the strong annihilation of the ip' into light hadrons is 
a two-step process: in the first step the ip' goes into two 
gluons in a 0++ or ^ state and an off-shell /ic(3526); 
in the second step the off-shell he annihilates into three 
gluons to produce light hadrons. Their argument implies: 
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(a) to leading order there is no strong decay amplitude 
for the processes ^' pn and tp' K*K] (b) the 12 % 
rule should hold for hadronic processes which take place 
via the QED amplitude only. As far as the second im- 
plication is concerned, the present data give different ra- 
tios between ^j]' and J/i/' decay for wtt*^ and tt+tt^ final 
states, both of which are electromagnetic processes. Here 
the event form factor effect is taken into account [s^] , the 
difference between two kinds of processes is still obvious. 
Besides the explanation for p7r puzzle, this model predicts 
a sizable i>' {TT~^Tr~ or 77) /ii(1170) branching fraction. 

In a recent paper [sl], Artoisenet, Gerard and Wey- 
ers (AGW) update and sharpen the above idea which 
leads to a somewhat unconventional point of view: all 
non-electromagnetic hadronic decays of the ip' goes via 
a transition amplitude which contain a cc pair. AGW 
provide two patterns for these two-step decays, the first 
is 

V-' ^ (2NP5) + (3g) . (4) 

The physics picture is as follows: the excited cc pair in 
the tp' does not annihilate directly. Instead, it spits out 
two non-perturbative gluons (2NPg) and survives in a 
lower cc configuration (1 or 1 ^) which then eventu- 
ally annihilate into 3g. The decays ip' {2Tr)J/tp and 
ijj' — > riJ/ip follow this pattern. The second pattern is 

^' ^ (3NP.9) + (25) , (5) 

where the lower cc configuration (0 or 0"*"^) annihilates 
into 2g. The only on-shell channel for this type of decays 
is Ip' — > (37r)77c, whose branching fraction is estimated as 
(1 — 2)% level. Anyway, the recent measurement from 
CLEOc [l^l provides the upper limit which is one order 
of magnitude below this theoretical prediction. Further- 
more, the substitution of one photon for one gluon in 
Eqs. (gl) and (O allows 

V^' ^ (2NP.9) + {2g) + 7 . (6) 

This pattern corresponds to on-shell radiative decays 
such as -0' ~> (7r+7r^)ryc7 and tp' —>■ rirjc^^ which could 
be larger than the observed ip' — > ?7c7 mode. 

Besides the above predications, AGW also estimate 

Bitp'^biT]) = (1.3 ± 0.3) X 10-3 , (7) 
i3(V'' ^ /ii7r°) = (1.9 ±0.4) X 10^3 , (8) 
BiJ/tP ^ hrj) « B{i;' 6177) « 1% . (9) 

All these are to be tested by experiments. 

6. Nonvalence Component Explanation 

Since the %p' is a highly excited state and close to the 
DD threshold, it is suggested [s^ that unlike the J/ip, 
the tp' may be an admixture of large nonvalence com- 
ponents in its wave function. The so-called nonvalence 



component indicates the additional gluon or the light 
quark-antiquark pair (or as in Ref. [5J], the ccg compo- 
nent and a quasi- molecular DD state), which makes tp' 
decays prominently distinctive from those of J/tp. The 
nonvalence component of the J/tp is expected to be less 
significant because it is the lowest state. In a following 
paper fi^, Chernyak uses this picture to explain the pir 
puzzle. He suggested that the valence and nonvalence 
strong contributions interfere destructively in p-n chan- 
nel and consequently cancel to a large extent in the total 
tjj' — > piT strong amplitude, while the role of nonvalence 
contributions is much less significant in J/tp — > pir. From 
this viewpoint, there is no deep reason for the experimen- 
tally observed very strong suppression of tp' pn, this 
is the result of a casual cancellation. 

Chernyak also tries to use the above idea to explain 
qualitatively other decay modes, such as VT, AP, PP, 
VV and BB decay. However, such an idea remains a 
pure speculation, and no concrete calculations have been 
performed up to now. 

7. S-D Wave Mixing Scheme 

The tp" is viewed as a Z^-wave charmomium state. Al- 
though it is primarily cc( I'^Di), its leptonic width in- 
dicates a contribution from mixing with S'-wave states, 
mainly the nearby tp{2^Si). This leads Rosner to propose 
that the small pn branching fraction in tp' decay is due 
to the cancellation of the 25 and ID matrix elements. 
By virtue of his scheme 

(pn\tp') = (pt:\2^ Si) cose -{pTr\l^Di) sine , 
{pTT\tP") = (p7r|23S'i) sine* + (pTrll^Di) cose , ^'^^^ 

where is the mixing angle between pure tp{2^Si) and 
tp{l^Di) states [11] and is fitted from the leptonic widths 
of the tp" and the tp' to be (12 ± 2)° [s^ , which is consis- 
tent with the coupled channel estimates [13, [11| and with 
the ratio of tp' and tp" partial widths to J/tpn~^Tr~ [soj . 
If the mixing and coupling of the tp' and tp" lead to com- 
plete cancellation of tp' — > pir decay {{pTT\tp') — 0), the 
missing pir decay mode of the tp' shows up instead as de- 
cay mode of the tp" , enhanced b y th e factor 1/ sin^ 9, the 
concrete estimation shows that [561 ] 

= (4.1 ±1.4) X 10-^ . (11) 

To test this scenario by the data collected at the tp" in 
e'^e- experiments, it has been pointed out [13, [Ml that 
the continuum contribution must be considered carefully. 
Specifically speaking, by Rosner's estimation, the Born 
order cross section for tp" — > pn is 

a|,?™^ = (4.8±1.9)pb, (12) 

which is comparable in magnitude to that of the contin- 
uum cross section, viz. 

af/™ -4.4pb. (13) 
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So what is observed is the coherent sum of the two am- 
phtudes. Based on the analysis of experimental data, 
it has been suggested that there be a universal phase 
between strong and electromagnetic amplitudes in char- 
monium decays. With this assumption, the strong decay 
amplitude interferes with the continuum amplitude ei- 
ther maximum destructively, e.g. for pn, Lurj and K*~^K~ 
or maximum constructively, e.g. for K*'^K^. The de- 
structive interference leads to the phenomena that the 
measured cross section on top of the resonance could be 
smaller than that off the resonance. Recent experimental 
results [iBIlll on pir, ujt] and K*'^K~ have demonstrated 
such interference pattern. This provides support to the 
prediction of Eg. (fTT|) . However, to uniquely determine 
B^ii^pTT, current available experimental data are insuf- 
ficient, the ■0" resonance must be scanned |62|]. So the 
quantitative test of Rosner's scenario remains to be the 
task of the future experiments on r-Charm factories. 

In the consequent study Wang, Yuan and Mo 

(WYM) extend the S-D wave-mixing scenario to PP fi- 
nal state, and give a constraint for -0" KgKl decay, 

0.12 ± 0.07 < 10^ X B{ij" ^ K"sKI) < 3.8 ± 1.1 , (14) 

which is compatible with both the BESII result B{ip" — > 
K'^Kl) < 2.1 X 10-4 at 90% C.L. fl^ and the CLEOc re- 
sult Biij" K^Kl) < 1.17 X 10-5 at 90% C.L^M. Ex- 
trapolating fore-study to all charmless decays [6J] , WYM 
found that this scenario could accommodate large non- 
DD decay of the ip" . Although the recent experimental 
measurement from CLEOc favors the zero non-DD cross 
section at tp" [6^ , the larger error could not rule out the 
existence of non-DD branching fraction at 10% level. 



C. Other Explanations 

Besides the models in two previous sections, more spec- 
ulations involving the pir puzzle are collected in this sec- 
tion. 



1. Final State Interaction Scheme 

Li, Bugg and Zou [6^ (LBZ) pointed out that the fi- 
nal state interactions (FSI) in J/ip and ip' decays give 
rise to effects which are of the same order as the tree 
level amplitudes, they argued that Jf-ip —> pi: is strongly 
enhanced by the a2p loop diagram, while the direct tree- 
production for pTT may be suppressed by the HHC mech- 
anism [2^ . The contribution of the aip loop diagram is 
much smaller than that of a2p loop for the J/V' — > Ptt, but 
they have similar strength for the ip' pir and may can- 
cel each other by interfering destructively. The similar 
apparent suppression for tp' K*K and /2a-' may also 
be explained by the K*K2 1 and hiir loop, respectively. 
Therefore, LBZ expected that FSI may provide a coher- 
ent explanation for all the observed suppressed modes of 



ip' decays. However, as remarked in Ref. ^7|i this inter- 
ference model does appear to have more assumptions that 
predictions. The only qualitative prediction by LBZ is to 
check whether the aip and K^K production rates are 
large for the ip' . So far, no such measurements have been 
reported. Nevertheless, useful information on aip and 
K^K could be obtained from the BES published data 
as shown in Refs. [13] and [l^. The lack of evidence 
within the invariant mass distribution plots (see Fig. 3 
and Fig.5 of Ref. [l3| or Ref. [li|) that the pi: recoiled 
against a p for events of tp' ~^ p^p^ir^ and that ir^K^ 
recoiled against a K*^ for events of ip' tt^tt^K^K^ 
suggests that they are unlikely to be the favored modes 
in -ip' decays. 



2. Large Phase Scheme 

Suzuki gave another explanation based on FSI for J /ip 
decays [38]. He performed a detailed amplitude analy- 
sis for J/ip decay to test whether or not the 
short-distance FSI dominates over the long-distance FSI 
in the J/tp decay. His result indicates that there is a 
large phase between three-gluon and one-photon ampli- 
tudes. Since the large phase cannot be produced with the 
perturbative QCD interaction, the source of it must be 
in the long-distance part of strong interaction, namely, 
rescattering among hadrons in their inelastic energy re- 
gion. Suzuki then performed the similar analysis for 
J/ip — > 0^0" decay, and obtained the same conclusion 
of the large phase [6^ . His analysis also shows that the 
exclusive decay rate at the J/ip is in line with that of 
the inclusive decay. This fact leads him to believe that 
the origin of the relative suppression of V'' — * 1^0^ to 
J /ip — > 1-0" is not in the J/ip but in the ip' . 

As to this large phase, in Ref. [KH, Gerard argued that 
this phase follows from the orthogonality of three-gluon 
and one-photon virtual processes. As a matter of fact, the 
conclusion of a large phase close to 90° between three- 
gluon and one-photon processes holds true for all two- 
body decays of J/ip, such as l+(r(90°)j6i, I'O" (106± 
10)° [Tdim 1"1~ (138^37)° [M,[nlJzilo-0" (89.6± 
9.9)° [61,1711, [73 and NN (89 ± 15)° ^U^. 

Nevertheless, as the ip' decays were analyzed, the ex- 
perimental data at first seemed to favor a phase close to 
180° [1^, in contrary to the expectation that the decay 
of J/ip and Ip' should not be much different. However, 
as pointed out by Wang et al. that the previous pub- 
lished data did not take the continuum one-photon pro- 
cess into account. Their reanalysis of the experimental 
data shows that the phase with value around —90° could 
fit iP' l-Q- data [m and ±90° could fit iP' O'O" 
data [t^ ]. The latter is confirmed by more recent re- 
sults by CLEO P]. Furthermore, this large phase also 
shows in the OZI suppressed decay modes of the ip" . 
In many decays modes of the ip" , the strong decay am- 
plitudes have comparable strength as the non-resonance 
continuum amplitude, the large phase around —90° leads 
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to destructive or constructive interference. Due to the 
destructive interference, the observed cross sections of 
some modes at the peak of the -0" are smaller than the 
cross section measured off-resonan ce |6Cl| . This is demon- 
strated by the data from CLEOc [2^ . 

If the large phase between the three-gluon and one- 
photon virtual processes is universal not only in J/?/; de- 
cays, but also in tp' as well as all charmonium or perhaps 
all quarkonium decays, then how this phase is related to 
the difference between the decays of J /ip and ^' remains 
to be an unanswered question. 

3. Mass Reduction Explanation 

In the study [zl] of radiative decays of 1 quarko- 
nium into 77 and 77', Ma presented a QCD-factorization 
approach, with which he obtained theoretical predictions 
in consistency with CLEOc measurement. The largest 
possible uncertainties in analysis are from the relativistic 
corrections for the value of the charm quark mass. Ma ar- 
gued that the effect of these uncertainties can be reduced 
by using quarkonium masses instead of using quark mass. 
As an example of such reduction, he provided a modified 
relation to the original 12% rule 

_ B{J/^ ^ pit) _ B{J/^J ^ e+e-) 

= (3.6 ±0.6)% . 

However, this value is much larger than the experimental 
result given in Table [H 

4-. Vector-meson-mixing Model 

Intending to give a comprehensive description of 
J/i/j two-body decays, Clavelli and Intemann (CI) pro- 
posed a vector-meson-mixing model in which the vector 
mesons {p, to ,(f> , J/ip) are regarded as being admix- 
ture of light-quark-antiquark state and charmed-quark- 
antiquark state [t^. The coupling of the J/^ to any 
state of light quarks is then related to the corresponding 
couphng of the p, w, and (j) to the same state. With few 
experiment inputs to determine the mixing parameters, 
CI calculated VP, PP, and BB decay rates for the J/ tp as 
a function of the pseudoscalar mixing angle. Most of the 
predictions agree with the experiment results at the order 
of magnitude level, but discrepancy is obvious for some 
channels, such as KgK^ final state @. CI also extended 
their model to the hadronic decays of the ijj' . Neverthe- 
less, their evaluations for B{J/ip ujtt'-') = 3 x 10~^ and 
B{'4'' wtt'^) = 3 X 10"'^ contradict with the present 
results (4.5 ± 0.5) x 10"'* and (2.1 ± 0.6) x IQ-^ [13, 
respectively. 

Starting from effective Lagrangian whereby nonet- 
symmetry breaking and pseudoscalar-meson mixing can 



be studied, Haber and Perrier parametrized the decay 
modes of J ftp — > PP (also for J/tp VV or rjc 
VP), J/tP VP (also for J/^P VT or rj^ ~^ VV), 
J/ij PPP (also for J/?/; ^ VVP or 77^ PPV), and 
?7c ^ PPP (also for J/tP PPV or 77^ VVP) [z3. 
The experimental data are used to determine these phe- 
nomenological parameters. In a later work, Seiden, 
Sadrozinski and Haber took the doubly Okubo-Zweig- 
lizuka suppression (DOZI) effect into consideration, and 
presented a more general parameterization of amplitudes 
for J/tp ^ PP decays [7g. With this form, one could 
easily derive the relative decay strength between different 
final states. However, it has also been noticed that un- 
der the most general circumstances, symmetry breaking 
introduces too many parameters for a general analysis to 
be useful. In order to reduce the number of new parame- 
ters considerably and make the analysis manageable, only 
special cases could be considered. 

TABLE II: Comparison of predictions 79] and measure- 
ments [l^l for the branching ratios (10~^) for tp' decays. The 
upper limits are presented at 90% C.L. 



VP 


Prediction Measurement 


pn 


1.3 


3.2 ± 1.2 




5.1 


10.9 ± 2.0 


K*+K- -\-c.c. 


1.3 


1 7+0.8 
-0.7 


UJTj 


1.2 


< 1.1 




6.3 


0+2.5 


<t>V 


1.6 


2 8+1-° 

^•°-0.8 




4.6 


3.1 ± 1.6 




2.1 


2.2 ±0.6 


pv' 


1.2 


1 9+^'' 
^■^-1.6 




3.8 


2.1 ± 0.6 


<^7r0 


0.01 


< 0.40 



A similar parameterization with mixing feature of the 
strong interaction mechanism was proposed by Feld- 
mann and KroU (FK) 79] for the hadron-helicity non- 
conserving J/tp and t/j' decays, but with a different inter- 
pretation from those put forth in Refs. [H, El, HE Izll ■ 
FK assumed that with a small probability, the charmo- 
nium possesses Fock components built from light quarks 
only. Through these Fock components the charmonium 
state decays by a soft mechanism which is modeled by 
J/ip-oj-cj) mixing and subsequent uj (or (jj) decays into the 
VP state. In absence of the leading-twist perturbative 
QCD contribution, the dominant mechanism is supple- 
mented by the electromagnetic decay contribution and 
DOZI violating contribution. FK argued that this mech- 
anism can probe the charmonium wave function at all 
quark-antiquark separations and feels the difference be- 
tween a IS" and a 25* radial wave functions. The node 
in the latter is supposed to lead to a strong suppression 
of the mechanism in ip' decays. With a few parameters 
adjusted to the experimental data, FK obtained a nu- 
merical description of the branching fractions for many 
VP decay modes of the J/4> and -0', which agree with the 
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measured branching fractions at the order of magnitude 
level, as shown in Table HIl Moreover, FK has extended 
their mixing approach to the rjc VV decays and ob- 
tained a reasonable description of the branching fractions 
for these decays while the rj'^ VV decays are expected 
to be strongly suppressed. 



D. Remarks 

From the brief retrospect of pn puzzle history, we 
could notice that in the early stage, theorists concen- 
trate on peculiar mechanism for special channels like pir 
and K*K, such as Hou and Soni's explanation. With 
the development of the theory and the progress of the 
experiment, theorists attempt to provide a more general 
scheme for charmonium decays, such as the work by Feld- 
mann and KroU. In fact, the charmonium decay is an 
interconnected system as a whole, a correct explanation 
of pn puzzle is expected to describe, quantitatively, or 
semi-quantitatively, the properties of all measured decay 
modes. 

At last, a few words about Fock state component. This 
concegt has been discussed in many works, such as by 
BK [4l|, CB dil, and FK [79|]. Even from the theoreti- 
cal explanation of experiments other than e'^e" collision, 
the non-perturbative Fock component is also indicated to 
exist. For example, the analysis of charmonia photopro- 
duction amplitude implies ^8Q] that the |cc) component 
is narrowly distributed in the transversal direction, while 
the cc pairs at larger separations may be part of higher 
Fock states which contain gluons and light quarks. 



III. EVALUATION OF RATIO OF i^' TO J/i/> 
DECAY 



In this section, we depict three approaches to estimate 
the ratio of ip' to J/ tp decay, one from theoretical analy- 
sis [H, the other two from experimental evaluation [sif . 



A. Theoretical Method 

Conventionally, the measured ratios of tp' to J/ip 
branching fractions for specific exclusive hadronic decays 
are compared with the naive prediction of pQCD, the so- 
called "12% rule". In the framework of pQCD Q, the 
■0 particles are considered to be nonrelativistic bound 
states of a charm quark and its antiquark. Their decays 
into light hadrons are believed to be dominated by the 
annihilation of the cc pair into three gluons. In order 
to annihilate, the c and c must have a separation of or- 
der l/rric, which is much smaller than the size of the 
charmonium state. Thus the annihilation amplitude for 
an S-wave state like J/ip and ip' must be proportional 
to the wave function at the origin, ip{r = 0) [l|. The 



width for the decay into any specific final state h consist- 
ing of light hadrons is therefore proportional to |'(/'(0)p. 
The width for the decay into e+e~ is also proportional 
to |-!/'(0)p. This leads to the simple prediction that the 
ratio of the branching fractions between tp' and J/V' is 
given by Eq. ([1]). 

However, this naive prediction suffers from several ap- 
parent approximations. Higher order corrections, which 
may not even be small, are not included in this calcula- 
tion. For example, a first order correction to the branch- 
ing fraction of J/'0 — > e~^e~ could be 50% of the lowest 
term if as{mj/^) ^ 0.2 (s^ is taken. The relativistic 
effect is also ignored. Since the mass difference between 
J/^jj and ^' is around 20% and < v^/c^ >^ 0.24 for 
J/f/;, this correction may be at the same level as the low- 
est order contribution |,82j . The inclusion of the finite 
size of the decay vertex will significantly reduce the ggg 
decay width of the J/tp [s^. Moreover, the effect of non- 
perturbative dynamics is neglected, the size of which is 
hard to estimate. Therefore, people may question the 
validity of the 12 % rule as a serious benchmark for com- 
paring experimental data. 



B. Experiment Estimation (I) 

The first experimental estimation is based on the as- 
sumption that the decays of J/tp and tp' in the lowest 
order of QCD are classified into hadronic decays (ggg), 
electromagnetic decays (7*), radiative decays into light 
hadrons {'ygg), and transition to lower mass charmo- 
nium states {ccX) [t^ . [sij . Thus, using the relation 
^{999) + ^{199) + '^(7*) + B{ccX) = 1, one can derive 
^{999) + ^{199) by subtracting B{'^*) and B{ccX) from 
unity. 

The calculated values of B{"f*) and B{ccX), together 
with the values used to calculate them are summarized in 
Table Hill As regards to ip' , two final states 777(25*) and 
hc{\^Pi)+X with faint branching fractions are neglected 
in our calculation. By deducting the contributions B{'^*) 
and B{ccX), we find that B{J/ip -> ggg) + B{J/i} 
-fgg) = (73.3 ± 0.5)% and S(^' ^ ggg) + S(^' ^ jgg) = 
(18.9 ± 1.3)%, then the ratio of them is 



Bj-ip' 999 + 199) 
B{J/ip -> ggg + igg) 



(25.7 ±1.7)%. (15) 



The above estimation is consistent with the previous 
ones [g^, ^ij. The relation between the decay rates of 
ggg and ygg is readily calculated in pQCD to the first 
order as 1831 



r( J/V; ^ 75.g) 16 



T{JI^ -> ggg) 5 (^^{7x1^) 



fl-2.9^ 



Using as{mc) — 0.28, one can estimate the ratio to be 
0.062. A similar relation can be deduced for the ip' de- 
cays. So we obtain B{J/^ — > ggg) ~ (69.0 ± 0.5)% and 
Bitp' ggg) ~ (17.8 ± 1.2)%, while the "25.7% ratio" 
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TABLE III: Experimental data on the branching fractions for 
J /tj) and ?/;' decays through virtual photon and to lower mass 
charmonium states used in this analysis. Most of the data are 
taken from PDG [IJ, except for B{J/ip, i/;' ^ 7* ^ hadrons), 
which are calculated by the product 7? • B{J/ip^ip' — > )j,'^fi~), 
with R = 2.28±0.04 [8J|. In estimating the errors of the sums, 
the correlations between the channels are considered [831 . 



and 



Channel 




'^\'r J 


7* — > hadrons 


(13.50±0.30)% 


(1.66±0.18)% 




(5.94±0.06)% 


(7.35±0.18)xl0-^ 




(5.93±0.06)% 


(7.3 ±0.8) X 10"^ 


r+T- 




(2.8 ±0.7) X 10"^ 


■y' ^ X 


(25.37±0.35)% 


(3.41±0.27)% 




(1.3±0.4)% 


(2.6±0.4)xl0~^ 


V'^ ■K~ J / ijj 




(31.8 ±0.6)% 






(16.46±0.35)% 






(3.09 ±0.08)% 


TT°J/tP 




(1.26±0.13)xl0"^ 


7XcO 




(9.2±0.4)% 






(8.7±0.4)% 


7Xc2 




(8.1±0.4)% 


ccX 


(1.3±0.4)% 


(77.7±1.2)% 



in Eq. (jlSp stands well for both ggg and jgg. Although 
Qg is considerably enhanced relative to Qh in Eq. ([1]), it 
is fairly compatible with the ratios for the K~^K~ and 
KgKl decay modes between ip' and J/ip, which are 



Qk+k- 
Qk%kI 
Qk%kI 



(26.6 ±4.5)% (CLEO) 
(28.8 ± 3.7)% (BES) , 
(32.2 ±5.2)% (CLEO) 



according to the recent results from CLEO and BES d, 
H, 0| ■ The relation in Eq. p5)) was discussed in the lit- 
erature as the hadronic excess in 'ijj' decays 0, l8l| ■ It 
implicates that while some modes are suppressed in ij)' 
decays, the dominant part of ^' through ggg decays is 
enhanced relative to the 12% rule prediction in the light 
of J/t/) decays. 



C. Experiment Estimation (II) 

The second approach for estimating Qh is to use 
the data on branching fractions for hadronic decays 
in final states containing pions, kaons, and protons 
that have already been measured for both the J /ip 
and the ij/ . They are tt+tt", K'^K^, pp, Tr+Tr^Tr*^, 
ppTT^, 2{TT+n-), 3(7r+7r-), 2(7r+7r-)7r°, 3(7r+7r-)7r°, 
2(7r+7r-7r"), 2{K+K-), K+R-Tr+n-, K+ R- tt+ tt- , 
K'^ K^2{tt^t:^), tt^tt^pp, ppn^n^TT^ and so forth. Us- 
ing the data compiled in Table HVl we have 
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^S(V''-^/0 = (1-70 ±0.23)% 



It follows that 



Qs = f: B{^' ->/.)/£ B{J/^p ^ /,;) 

j=l / i=l 

= (11.8 ± 1.7)% . 



(16) 



TABLE IV: Branching fractions for ip' and J/ip decays, and 
Qh values are also calculated. For ip' — > tt+tt", K~ de- 
cays, the results are the weighted average of measurements 
from CLEO H and BES [§1], while for J/V- -k+ti'-k^, the 
result is the weighted average of measurements from BES 0] 
and BABAR 0|. Except for values with f, all others from 

PDG na. 



Final state 




B*(25)(10-") 


Qh (%) 




0.147 ±0.023 


0.08 ± 0.05t 


5.6 ±3.5 


K+K- 


0.237 ±0.031 


0.63 ± 0.07t 


26.5 ±4.5 


PP 


2.17 ±0.08 


2.65 ±0.22 


12.2 ± 1.1 




21.2 ± l.Olf 


1.68 ±0.26 


0.79 ±0.13 




1.09 ±0.09 


1.33 ±0.17 


12.2 ± 1.9 


2(7r+7r~) 


3.55 ±0.23 


2.4 ±0.6 


6.8 ± 1.8 


3(7r+7r-) 


4.3 ±0.4 


3.5 ±2.0 


8.1 ±4.7 


3(7r+7r-)7r0 


29 ±6 


35± 16 


12.1 ±6.1 


2(7r+7r-)7r'' 


33.7 ±2.6 


26.6 ±2.9 


7.9 ± 1.1 


2(7r+7r-7r'') 


16.2 ± 2.1 


45 ± 14 


27.8 ±9.4 


2{K+K') 


0.78 ±0.14 


0.60 ±0.14 


7.7 ±2.3 




7.2 ±2.3 


7.2 ±0.5 


10.0 ±3.3 




12.0 ±3.0 


12.4 ± 1.0 


10.3 ± 2.7 


K+K-2{TT+-n-) 


4.7 ±0.7 


18 ±9 


38.3 ± 20.0 




6.0 ±0.5 


6.0 ±0.4 


10.0 ± 1.1 




2.3 ±0.9 


7.3 ±0.7 


31.7 ± 12.8 



Y,l3iJ/iJ ~* fr) = (14.46 ±0. 



Some remarks are in order here. First, we know that 
most of the multihadron final states in fact include sums 
of several two-body intermediate states, so the Qs is 
not the exact ratio of i/j' to J/ijj inclusive hadronic de- 
cay rates, but the ratio on average of the exclusive de- 
cay channels as measured to date. In another word, 
Qs represents a mixed effect which may deviate notice- 
ably from the component Q-values. For example, the 
decay ip' — > 7t'^tt~ K~ can proceed through inter- 
mediate state K*{892fKl{U30)° + c.c, whose Q = 
(2.9 ± 1.3)%, is greatly suppressed {l9\ comparing with 
Q — (10.0 ± 3.3)%. Second, we notice that the results 
obtained by two estimations vary considerably; further- 
more, by virtue of Table IIV| it is obvious that many 
Q- values deivate from Qs significantly while on average, 
Qs is similar to pQCD Qh- It would be an intriguing 
problem that with more and more data and higher accu- 
rate measurements if Qs and Qg could approximate with 
each other and be consistent with Qh, or if Qs and Qg 
still deviate from each other prominently. 
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IV. COMMENT ON THE pvr PUZZLE 
A. Theory Aspect 

For all explanations involving pir puzzle, there is a ba- 
sic assumption that the non-relativistic potential model 
delineates the physics of charmonium decays to a good 
approximation. However, this assumption indeed re- 
quires examination in detail. First, as we have noticed 
in subsection IIII Al several corrections, should be added 
for the decay ratio between -0' and J/i^- 

Second, we see the effect of non-perturbative dynam- 
ics is neglected, which is crucial for charmonium decays. 
Actually from various explanations of pir puzzle, or more 
generally from the phenomcnological explanation involv- 
ing pQCD, we note that certain non-perturbative effect 
or nonlinear effect must be incorporated one way or the 
other in order to recount the experimental data. But this 
kind of effect could hardly been included in the present 
non-relativistic potential model. 

Third, if the S-D wave mixing scenario holds as the 
solution of pTT puzzle, then the matrix element of D wave 
to light hadrons would be very large, which can hardly be 
accommodated in the potential model, or any other cur- 
rent theory. It indicates that the current understanding 
of charmonium decays may not be complete. 

We envisage that new development of the theory 
should take into account the following features involving 
charmonium decays 

• the mass effect for different charmonium states, 
such as J/ip, ip' , ip" and so on; 

• the non-perturbative or non-linear effect on the res- 
onance decays; 

• reasonable description for the known features of 
charmonium spectroscopy; 

• quantitative consistence (with reasonable high ac- 
curacy) with the existing experimental measure- 
ments. 

Herein it is also important to distinguish the quarko- 
nium states (theoretical states) and mass eigenstates 
(physical states) [sil . [s^ . From a fairly theoretical point 
of view, if both states, say, J/ip and tp' did dominantly 
annihilate into three gluons, they would mix and could 
thus not be the putative quarkonium states. In a non- 
relativistic potential model, for example, the ip' is simply 
a radial excitation of the J/4'- This is a well defined pic- 
ture in which J/tp and ip' are orthogonal states. If the 
annihilation into three gluons could be treated as a "per- 
turbation" to the non-relativistic potential, then clearly 
the unperturbed states would mix and rearrange them- 
selves into orthogonal mass eigenstates. The QCD dy- 
namics may be such that the physical states, presumably 
mixtures of the theoretical cc states {n?^~^^Lj), are so 
built up that one of them strongly annihilates into three 



perturbative gluons while the other does not. Mixing of 
the l^Si and 2^Si states via three perturbative gluons 
has little effect on the charmonium mass spectrum, but 
may be crucial for the decay pattern. 



B. Phenomenology Aspect 

One may remember that at the early stage of the dis- 
covery of a narrow state J/ip, the cc system was hailed as 
the Hydrogen atom of Quantum Chromodynamics, with 
the implied hope that the study of the newly discov- 
ered system could shed as much light on the dynamics of 
quark-antiquark interactions as the study of the Hydro- 
gen atom had on Quantum Electrodynamics. But one 
may also notice the historical fact, even before Bohr's 
theory, Balmer series had been discovered for long, and 
the famous Rydberg formula had also been proposed, 
which laid a solid foundation for further theoretical im- 
provement. If we are conscious of the more complicat- 
edness of charmonia system comparing with Hydrogen 
atom, we may prepare for more hard and meticulous 
works. As a first step, it is necessary to develop a reli- 
able and extensively applicable phenomcnological model 
(PM). 

The advantage of PM lies in that a well-defined 
PM contains few experimentally determined parameters 
which have clear physical meaning; moreover, with only 
few parameters determined from experiment, PM could 
produce concrete results which can be directly confirmed 
or falsified by experiment and may guid experimental 
searches. Such a model has a good relation with elemen- 
tary principle of the theory, and if correct, can be used 
for further theoretical refinement. This point is note- 
worthy for the time being, since the general QCD can 
hardly provide solutions for special problems; conversely, 
we have to establish certain effective empirical model to 
advance our understanding for generic QCD principle. 

Here we would like to mention few ideas of PM, which 
have or intend to provide a general description for char- 
monium decay. 



1. Mixing Model 

As we have noticed in subsection III C 41 whatever CI 
model or FK model, they could yield definite predic- 
tions for experimental test, and therefore provide clues 
for further development. In addition, according to Haber 
and Perrier's parameterization scheme [T^Ij could get 
the decay rate relations between different channels, and 
the proportions of different interaction amplitudes, all of 
which are useful information to understand the dynamics 
of charmonium decays. 

The s pec ial feature of Rosner's S- and D-wave mixing 
scheme [56| is that it is simple, and it works both for 
suppressed and enhanced decay modes, and moreover, it 
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connects the J/^p, ip' and tp" decays together, and give 
predictions with httle uncertainty. 



2. Effective FSI 

Miller once discussed the connection between the 
strong-coupling approximation to quantum chromody- 
namics and nuclear properties observed at low and 
medium energies and momentum transfer [s^]- He sug- 
gested that the strong-coupling (corresponds to long- 
range interaction) version of QCD does reproduce the 
salient feature of the meson-baryon picture of low mo- 
mentum transfer in nuclear physics. The derivation of 
nuclear physics indicates that quark aspects of ordinary 
nuclei are hidden in the hadronic degrees of freedom. Re- 
versely, one may imagine that a quasi-meson or quasi- 
baryon structure could be formed in non-perturbative 
hadronization process. The meson-baryon picture in nu- 
clear physics could be utilized, with some modification, as 
an effective FSI theory, or molecular-model theory. The 
recent molecular interpretations [H, [s^ for the newly 
found state y(4260) could be treated as such an effective 
FSI theory. 

Furthermore, the quasi-particles could even be real 
particles, and then with the residual strong force between 
the quarks inside the quasi-particles, some multiplets can 
be formed as the quasi-particle and anti-quasi-particle 
bound states. This idea was first put forth to explain a 
lots of low mass baryon-antibayron enhancements newly 
found by CLEO, Belle and BaBar collaborations [9Q] . 



3. Glueball and Hybrid 

We have seen that the J/i/j-glueball mixing scheme is 
the first explanation proposed to explain the ott puzzle. 
With the implication from lattice calculation [91| , Suzuki 
once proposed a glueball-'0' mixing scheme to explain the 
excess hadronic decays at the ip' jl2] ■ In fact, gluodynam- 
ics is always a tantalizing domain for theorists. Recently, 
some lattice evaluations suggest [o^ [oj] that the masses 
of some hybrid states could be low enough to be degen- 
erated with charmonium states, such as J/ip and/or XcJ- 
Unlike Glueball, these ccg hybrid states could be very 
narrow, so it may be very stimulating and noteworthy to 
search for and confirm these kinds of states. 



C. Experiment Aspect 

Physics is a science of experiment. Physical facts are 
bases for theoretical development and also criteria for 
checking theoretical hypotheses. It is a prominent fact 
that the early analyses on pn puzzle based on meager ex- 
perimental data often lead to unsatisfactory, sometimes 
premature, inferences, which were washed out easily by 



later accurate data. In fact, the current knowledge con- 
cerning the ip' decays from experiments is still rather lim- 
ited, even summing all charmless channels presented in 
PDG [27II, the total branching fraction is less 2%. Such 
situation prevents us from laying down a solid founda- 
tion for elementary dynamics exploration. However, the 
estimation of Q value, as discussed in Section [1111 affords 
us some clues concerning the exploration of charmonium 
decay dynamics. Since many suppressed channels have 
been found, especially those such as p7r which is greatly 
suppressed in ip' decay, and if the Qh really represents 
the averaged value of inclusive hadronic decay, the esti- 
mation of Qg indicates that either lots of enhanced de- 
cays are not discovered, or some particular decays only 
present in ip' , or both cases exist. Therefore, systematic 
experimental study of ■0' decays is anxiously awaited. 

Moreover, study should be carried out not only on the 
ratios of ip' to J /ip decays, but also on the other ratios 
such as those between 77^ and 77c [ZQ], and those between 
i>" and J/V'[53|) and/or many other ratios between differ- 
ent resonances for the same channel or between different 
channels about the same resonance [sj] . All these studies 
will shed light on the understanding of charmonium de- 
cays. From our point of view, the progress would be more 
likely obtained from the analysis of the new experimen- 
tal results, rather than from the inspiration of general 
theoretical principle. 



V. SUMMARY 

In this paper, we present a general review on the study 
of the pTT puzzle. Although there is still no satisfactory 
explanation for all existing experiment results, some en- 
lightenment ideas have been put forth. In addition, we 
also discuss three methods of estimating the ratio of the 
branching fractions in J /ip and ip' decays. In the light of 
the present theoretical and experimental status, we argue 
that it is important to explore the potential models from 
a new point of view, it is necessary to search for and/or 
construct an effective phenomenological model, and it is 
especially crucial to perform systematic measurements of 
various charmonium decays. 

As we know physics in the charm energy region is in 
the boundary domain between perturbative and nonper- 
turbative QCD. Recently the observed hadronic decays of 
charmonium may give new challenges to the present theo- 
retical understanding of the decay mechanisms. Many of 
the new charmonium states observed by Belle and BaBar 
which can hardly be accommodated by potential models 
have led to new theoretical speculations about charmo- 
nium spectroscopy and novel production mechanism [osj . 

Experimentally, the expected large data sample from 
CLEOc in the near future, and even larger sample from 
forthcoming detector BESIII will open for us a new era 
for charmonium dynamics study, even though we may not 
obtain a throughly revolutionary theory, we could acquire 
more profound understanding for the existing theory, at 



12 



the same time we could expect a brand-new comprehen- 
sion for the pTT puzzle. 
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